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ABSTRACT 

The Seyfert galaxy NGC 1097 has an extended neutral hydrogen disk, a companion, a prominent 
bar and a luminous circumnuclear starburst "ring" . Magnetic fields as revealed by nonthermal radio- 
continuum emissions correlate well with the optical barred spiral structure on large scales, have a gross 
enhancement overlapping with the optical/infrared "ring", and show a trailing swirl around and within 
the "ring" . We propose a scenario of bar-excited long-trailing fast magnetohydrodynamic (MHD) density 
waves at the modified inner Lindblad resonance (mILR) , physically identified with the outer rim of the 
"ring" . These sustained outgoing long- waves are bounced back by the Qj\/-barrier in the form of incoming 
short-trailing waves. The damping of these waves deposits a negative angular momentum into the 
magnetized circumnuclear gas disk. Thus, gas materials spiral inward, bring in frozen-in magnetic flux, 
and accumulate inside the mILR to create a circular zone of high density and magnetic flux vulnerable 
to massive star formation. Depending on the wave damping efficiency, this process may simultaneously 
sustain a net mass inflow across the "ring" and toward the nucleus. A wavelet analysis on a Hubble 
Space Telescope image of central NGC 1097 shows a distinct two-arm spiral structure extended down to 
the nucleus as a strong evidence for circumnuclear MHD density waves. We predict that magnetic-field 
observations with improved sensitivity and resolution would reveal a specific correspondence between 
circumnuclear optical and magnetic field spirals much as those known to exist on large scales in nearby 
spiral galaxies, including NGC 1097. 

Subject headings: accretion, accretion disks — magnetohydrodynamics — galaxies: NGC 1097 — 
galaxies: spiral — galaxies: starburst — ISM: clouds 



1. INTRODUCTION 

In barred spiral galaxies, luminous circumnuclear rings 
of star bursts (Buta 1986; Telesco 1988; Barth et al. 
1995) commonly appear. Some galaxies with active nu- 
. clei (AGNs) are also known to contain circumnuclear star- 
burst rings (Simkin et al. 1980; Arsenault 1989; Barth 
et al. 1995). The origin of such circumnuclear activities, 
monitored in a wide spectrum including optical, infrared 
and radio-continuum bands etc, remains to be understood. 
(Lynden-Bell 1969; Shlosman et al. 1990). Even in the 
absence of self-gravity and magnetic field, numerical sim- 
ulations (Piner et al. 1995) have produced gross features 
of circumnuclear rings in a fixed bar-potential. While the 
driven flux of angular momentum is independent of the de- 
tails of the disk modeling, the exclusion of density waves 
that involve self-gravity will however affect the removal 
efficiency of disk angular momentum and the absence of 
magnetic field would at least make the formation of mas- 
sive stars in the circumnuclear "ring" unlikely (Shu et al. 
1987; Elmegreen 1994). In this Letter, we advance a physi- 
cal scenario for bar- or satellite-excited fast magnetohydro- 
dynamic (MHD) density waves (FMDW) in a magnetized 
self-gravitating disk (Fan & Lou 1996; Lou & Fan 1998; 
Lou et al. 2001) at the modified inner Lindblad resonance 
(mILR) to illustrate the processes of spiral FMDW prop- 



agation, angular momentum transfer, forming circumnu- 
clear rings, and wave-induced net mass accretions across 
the ring and toward the nucleus, and to comprehend op- 
tical, CO, Hi, infrared and radio-continuum observations 
of the barred spiral galaxy NGC 1097 in dynamic, mor- 
phologic, and diagnostic contexts (Rickard 1975; Meaburn 
et al. 1981; Ondrechen & van der Hulst 1983; Hummel et 
al. 1987; Gerin et al. 1988; Ondrechen et al. 1989; Barth 
et al. 1995; Beck et al. 1999; Kotilainen et al. 2000; 
Emsellem et al. 2001). While we focus on NGC 1097, 
the scenario is generally applicable to circumnuclear spi- 
rals and starburst rings in magnetized barred spiral galax- 
ies (e.g., NGC 6951, NGC 2997 etc.). This conceptual 
framework also bears broad implications to accretions and 
MHD processes in other astrophysical disk systems such 
as a protostellar disk system. 

With a high stellar velocity dispersion ( > 100 km s _1 ) 
and a low gas temperature, the central circumnuclear re- 
gion of a spiral galaxy contains a thin stable 1 gas disk 
within a few kiloparsecs (kpc). The gas rotation curve in 
an extended radial range is determined by the total mass 
distribution in the galaxy (Meaburn et al. 1981; Blackman 
1981; Gerin et al. 1988; Ondrechen et al. 1989) and the 
background magnetic field in the circumnuclear disk plane 
of a strength > several tens of [iG (Beck et al. 1999) 



1 Stability against bar-type instabilities is meant. Starbursts involve gravitational collapses on subscales along an accumulated circular zone. 

2 It is an idealization to posit a purely azimuthal background magnetic field, based on which we perform analysis for spiral MHD density waves. 
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is taken to be axisymmetric and azimuthal 2 to avoid the 
magnetic-field winding dilemma (Roberts & Yuan 1970). 
This circumnuclear magnetized gas disk is embedded with 
molecular clouds (Gerin ct al. 1988) with a mean-free path 
l c and a velocity dispersion v c . When disturbed, spiral 
FMDWs involve collectively the magnetized gas disk cm- 
bedded with clouds. The large-scale barred spiral structure 
outside the central gas disk rotates at a pattern speed fl p 
and gives rise to a time-dependent external periodic gravi- 
tational potential <fi E as felt by the central magnetized gas 
disk. In cylindrical coordinates (r, 6, z), our theory is de- 
veloped using the standard two-dimensional equations for 
coplanar MHD density waves in the central magnetized 
gas sheet (Lou & Fan 1998) at z — with <j> E varying 
slowly in r as an inhomogeneous driving term (Goldrcich 
& Tremaine 1978a,b; Yuan & Cheng 1991). 

2. KEY RESULTS OF THEORETICAL ANALYSES 

Without such cj> E driving and away from the resonances, 
the dispersion relation for free FMDWs in the WKBJ or 
tight-winding regime (Shu 1970; Fan & Lou 1996) is 
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where v 2 mimics 
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re" due to random 
cloud velocities, u is the angular frequency in an in- 
ertial frame of reference, m > an integer for the 
number of spiral arms, the disk angular speed, k = 
[(2Q/r)d(r 2 Q) / dr] 1 / 2 the epicyclic frequency, G the gravi- 
tational constant, no the background surface mass density, 
k the radial wavenumber, Cs the sound speed, and Ca 
the Alfven speed (Lou & Fan 1998). The corresponding 
FMDW amplitude equation (Fan & Lou 1999) is 
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where v r (r) is the magnitude of the radial velocity v r . 
Equation (1) is quadratic in |fc| and contains the familiar 
short- and long-branches of FMDWs (Lou & Fan 1998). 
For FMDWs, it follows that Bgfi = f-obg, where \x is the 
surface mass density perturbation, and bg and Bg are the 
perturbation and background azimuthal magnetic fields; 
the enhancements of \x and bg are therefore in phase. 

The stability against axisymmetric ring fragmentation 
requires the MHD version of Toomre's Q parameter Qm = 
k(C a + Cg + v 2 ) 1 / 2 /(nGfio) > 1 in a magnetized self- 
gravitating gas disk (Lou & Fan 1998). The radial group 
velocity of FMDWs is 

Ouj_ „ _ sgn(fc)[(C72 + C 2 + v 2 c )\k\ - ttGH 
dk 
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with Fq > and Fq < for outward and inward propaga- 
tions. By our sign convention, k > and k < represent 
leading and trailing spirals. Inside corotation of FMDWs, 
a wave packet travels away from the disk center for long- 
trailing and short-leading FMDWs; a wave packet travels 
toward the disk center for long-leading and short-trailing 
FMDWs. Outside corotation of FMDWs, the moving di- 
rections of wave packet reverse for all these FMDW types. 
The angular momentum flux carried by a FMDW is 
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with the surface densities of angular momentum J F = 
m/j,o\v r \ 2 /[2(cu — mfl)], of energy £ F = uuJ F /m and of 
wave action N F = J F /m. Inside and outside corota- 
tion, J F is negative and positive, respectively. By (2) 
and (4), the angular momentum flux of free FMDWs is 
conserved. The corotation at u — mCl(r c ) — is forbid- 
den to wave access by the QM-barrier when Qm > 1, and 
in the WKBJ regime, the slightly modified Lindblad reso- 
nances occur at T = k 2 - (u> - mfl) 2 + m 2 C 2 A f(r)/r 2 = 
with f(r) ~ O(l) being a dimcnsionlcss analytic expres- 
sion. Outside the QM-barrier, long-FMDWs exist only 
within the two modified Lindblad resonances, while short- 
FMDWs propagate within and outside the two modified 
Lindblad resonances as in the hydrodynamic case (Gol- 
dreich & Tremaine 1978a). 

By a massive bar, a satellite, or a large-scale spiral struc- 
ture, an external potential <fi E (r, 8, t) felt at the central gas 
disk excites FMDWs at the modified Lindblad resonances 
rjf. Around tm, one may write x = (r — rM)/rM and 
r = Qm x with sgn Qm = ±1 for modified inner and outer 
Lindblad resonances (mILR and mOLR). In the WKBJ 
regime, the set of inhomogeneous FMDW equations for 
small x may be reduced to the familiar form in terms of 
the disk self-gravity potential <f) as 
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where 
Pm = [ 

2mil (f> E / (mfl — u>). The solutions of (5) involving Airy 
functions with proper energy flow directions are known; 
the physics of these solutions studied previously (without 
magnetic field) in contexts of exciting large-scale galac- 
tic density waves (Goldreich & Tremaine 1978b; Yuan & 
Kuo 1997) and of planetary rings (Goldreich & Tremaine 
1978a, 1982) may be extended and applied to a circum- 
nuclear magnetized gas disk. By (5), the magnitude of a 
FMDW is proportional to the strength of * that involves 



(5), long-trailing FMDWs are domi- 
(f> E at mILR and mOLR; these waves 



In the circumnuclear domain under consideration, observed magnetic field deviations are significant 
only indicative. The effectiveness of this approach may be questionable in the fully nonlinear realm. 



By solutions of 
nantly excited by 
exist between the two resonances and propagate toward 
corotation after their generations at mILR and mOLR. 
Reflections at the Q M-barrier transform these long-trailing 
FMDWs to short-trailing FMDWs that travel away from 
corotation and continue across the mILR and mOLR if 
dissipations were not strong. In the case of mILR, long- 
trailing FMDWs sustained by cf> E carries a negative an- 
gular momentum toward corotation; a subsequent reflec- 
tion at the QM-barrier gives rise to short-trailing FMDWs 
again carrying a negative angular momentum that may 
reach the disk center until being disrupted by small-scale 
( < a few pc) processes around the nucleus. 



the circumnuclear spiral arms and starburst 
"ring" in ngc 1097 

NGC 1097 is a barred spiral galaxy with a declination of 
—30° and an estimated inclination ranging from ~ 37° 

The results of a perturbation analysis is 
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to 57°. It is of Scyfcrt type with an AGN known as a 
LINER (Barth et al. 1995) and may host a supermassive 
black hole (SMBH) of ~ 2 x 1O 1O M (Rickard 1975). Es- 
timates of its distance range from 11.6 — 17 Mpc. It has 
a bright circumnulear starburst "ring" with a diameter of 
~ 18" and a far infrared (FIR) luminosity of <~ 3 x 1O 1O L 
(Tclesco 1988; Gerin et al. 1988), forming stars at a rate of 
~ 5M yr _1 (Hummel et al. 1987). The rotation curves 
were determined for r ~ 4" — 20", using H a 6563A and 
[Nil] 6584A (Meaburn et al. 1981; Blackman 1981) as well 
as CO (J=l - 0) (Gerin et al. 1988), and for r ~ 20" - 6', 
using neutral hydrogen Hi 21 cm despite a strong depletion 
in the bar region (Ondrcchcn et al. 1989). From r <~ 10" 
to 20", the inferred rotation speed V$ <~ 300 km s^ 1 re- 
mains roughly constant. The two dark dust lanes origi- 
nating around <~ 4" — 5" follow along the inner edges of 
the two spiral arms and extend almost continuously out 
to the bar (Rickard 1975) with large arm pitch angles of 
<~ 50° there. Shown in Fig. 1 is our wavelet analysis on an 
HST image of central NGC 1097 that reveals two distinct 
circumnuclear spiral arms extended inward to < 1.5". 

Synchrotron radio-continuum emissions (Ondrechen & 
van der Hulst 1983; Hummel et al. 1987) at 20 cm reveal 
that the intensity contour closely follows the dust lanes 
along the bar and merges into the central circular region 
around r ~ 30" to 35". The total radio intensity ring over- 
laps with the circumnuclear optical/infrared "ring" (Hum- 
mel et al. 1987; Telesco 1988), implying strong magnetic 
fields and profuse relativistic cosmic-ray electrons. Polar- 
ized radio emission at 22, 18, 6.2 and 3.5 cm (Beck et al. 
1999) indicates that regular magnetic fields are parallel 
to the two dust lanes along the bar and gradually swirl 
into a trailing pattern in the central region. Relative to 
the circumnuclear "ring" , magnetic fields are inclined by 
~ 50°. Despite the limited resolution, it is already in- 
dicative that a trailing spiral magnetic-field pattern inside 
the "ring" continues toward the nucleus. By this evidence 
and our Fig. 1, it is almost certain that spiral FMDWs 
(Fan & Lou 1996; Lou & Fan 1998) play a significant dy- 
namic role in the circumnuclear region of a barred spiral 
galaxy, and magnetic fields provide a valuable diagnostics. 
The fact that in the "ring" of NGC 1097, the radial and 
azimuthal components of magnetic field are comparable 
(Beck et al. 1999) means that the magnetic-field swirl pat- 
tern is not so tightly wound there. This perhaps relates 
to long-trailing FMDWs that are excited and sustained 
by 4> E at the mILR, and spiral FMDWs involve a nonax- 
isymmetric radial magnetic field b r . Fig. la of Beck et al. 
(1999) with 15" resolution shows a central swirl with some 
axisymmetry, while their Fig. 2 with 6" resolution shows 
a central spiral pattern with nonaxisymmetry. The ap- 
parent central "axisymmetric swirl" of their Fig. la most 
likely results from a larger beam size. In terms of orienta- 
tions along magnetic spiral arms, the tendency that mag- 
netic field of central NGC 1097 seems to be axisymmetric 
instead of bisymmetric within the circumnuclear "ring" 



region may be understood from the intuitive picture of 
distorted unaligned magnetic-field rings for FMDWs (Lou 
& Fan 1998) similar to that of distorted stellar orbits for 
density waves (Kalnajs 1973). 

For NGC 1097, the mean magnetic field Bg threads 
through the central thin Hi gas disk embedded with H 2 
molecular clouds (Gerin et al. 1988). Spiral FMDWs in- 
teract with the magnetized Hi gas 3 and H 2 clouds together. 
Random motions of H 2 clouds provide an "effective pres- 
sure" and a necessary mechanism for FMDW damping. 
The mILR is physically identified to be located at r ~ 10" 
somewhat outside the radius 4 of the circumnuclear "ring" . 
For to = 2, 17 at mILR is 2 + \/2 times Q p in a flat rota- 
tion curve, the corotation is located around 30" < r < 35", 
and the mOLR would be around 50" < r < 60". The lat- 
ter two estimates are only indicative, as both stellar and 
magnetized gas disks are involved at larger r. 

Our emphasis is on FMDW processes around the mILR. 
Long-trailing FMDWs are excited and sustained at the 
mILR by cf) E and propagate toward corotation. They 
are reflected by the QAf-barrier around corotation and 
travel back toward the mILR in the form of short-trailing 
FMDWs. Inside corotation, both long- and short-trailing 
FMDWs carry negative angular momenta (Fan & Lou 
1999). Owing to dissipation, these FMDWs damp and 
deposit in the disk their negative angular momentum as 
they travel. As the disk angular momentum is persistently 
reduced, disk material outside the mILR gradually spirals 
inward, bringing along magnetic flux meanwhile. As this 
process persists, disk materials (Hi and H 2 etc) and mag- 
netic flux would accumulate inside the mILR. Regions of 
high gas density and enhanced magnetic flux naturally fa- 
vor births of bright young massive stars as well as star 
clusters (Elmegreen 1994); this presumably gives rise to 
a circumnuclear starburst "ring" just inside the mILR. 
Outside the mILR, the strength of a long-trailing FMDW 
should be stronger than that of a returning short-trailing 
FMDW owing to wave damping. 

For order-of-magnitude estimates, we take v c ~ 
35 km s _1 (Gerin et al. 1988) 5 and l c ~ 100 pc. The 
gross turbulent dissipation coefficient v c would be v c <~ 
v c l c /3 = 3 x 10 26 cm 2 s _1 . For a magnetic field of 
B ~ 40/zG (Beck et al. 1999), a Hi gas column density 
of < 10 21 cm~ 2 (Ondrechen et al. 1989), a H 2 gas col- 
umn density ct ~ 6 x 10 22 cm~ 2 (Gerin et al. 1988) 
and a disk thickness h < 100 pc, the Alfven speed is 

Ca ~ 4 x 10 5 cm s _1 . For a gas temperature T < 100 K°, 
the sound speed is Cs < 10 5 cm s _1 . A rule of thumb 
for active star formation in a magnetized gas disk is 
27tG 1 / 2 E/.B > 1 with S being the surface mass density 
(Shu 2001, private communications). Given the above pa- 
rameters, 2ttG 1 / 2 'E/B ~ 8. The Qm estimate is uncertain, 
as the corotation region 30" < r < 35" for two-armed spi- 
ral FMDWs may involve both stellar and magnetized gas 
disks; we take Qm > 1- Right at the mILR (r <~ 10"), the 



3 A partial ionization and sufficiently frequent collisions between neutral Hi gas and charged particles are presumed here. 

4 For the Cassini division in the Saturn ring system (Goldrcich & Trcmainc 1978a), particles outside the Mimas' 2:1 ILR have been cleared 
by density waves over 5 billion years. In the disk region outside the circumnuclear starburst ring of a barred spiral galaxy, there may be some 
mechanisms in operation to replenish the gas disk to avoid an outright gap. 

5 There seems to be no escape for a high Mach number of molecular clouds in the gas disk. Resonance processes might continuously convert 
FMDW energy into random cloud motions. 
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short wavelength is \ s = {C\ + Cf + v*)/{G(j, ) ~ 300 
pc, while the long wavelength A^ would be infinite; at 
the reflection point of the QM-barrier (r ~ 30"), they be- 
come equal \ s = \ L = 2{C\ + C*| + vl)/(Gfi ) ~ 600 
pc. For a mean radial wave lengthscale of Am ~ 450 pc, 
the damping timescale would be r c ~ }? M jv c = 6 x 10 15 
s. The timescale r w for a FMDW to go from the mILR 
to QM-barrier and back to the mILR is estimated to be 
t w <~ 3 x 10 15 s. Thus, a sizable fraction of negative an- 
gular momentum carried by long-trailing FMDWs excited 
and sustained at the mILR would have been deposited in 
the gas disk as a weaker short-trailing FMDW returns to 
the mILR. As the disk angular momentum is reduced, gas 
materials with frozen-in magnetic flux from outside tend 
to accumulate just inside the mILR to form a circumnu- 
clear starburst "ring" . An estimate of net mass inflow rate 
M involves several uncertain aspects because the magni- 
tude of excited FMDWs relates to the strength of <j) E and 
the wave damping distribution depends on the effective 
viscosity as well as some nonlinear effects. By taking a 
\v r \ <~ 10 km s _1 and a \u> — mO|r <~ 30 km s _1 in (4), an 
upper limit of M may be estimated as M < IOMqvt -1 . 
For a stronger cf) E , nonlinear wave and damping effects 
may become important and this upper limit may be raised. 



superposed with the stronger long-trailing FMDW. If one 
were to regard the large magnetic-field pitch angle as a 
non-frozen-in steady state, which we think unlikely, then 
the upper limit of M would be ~ lM Q yr _1 (Beck et al. 
1999). The upper limit for M induced by turbulent viscos- 
ity with a flat rotation curve is of the order of ~ lM Q yr _1 . 
Some simulations also show that viscous torques are neg- 
ligible with respect to gravitational torques several kpc 
away from the nucleus, when a relatively strong bar or 
tidal interaction is present (Combes 1991). In terms of 
feeding active star formation, the relative importance of 
viscous and gravitational torques around circumnuclear 
rings remain unresolved. 

In our scenario, the effectiveness of wave damping holds 
the key for different outcomes. For a weak damping, 
FMDWs excited at the mILR eventually travel to the 
nucleus. For a strong wave damping, disk angular mo- 
mentum is reduced primarily outside the mILR. For an 
intermediate wave damping, in addition to a circumnu- 
clear "ring" formation at the mILR, remnant short-trailing 
FMDWs also travel across the mILR and be partially en- 
hanced by the geometric converging effect (Montenegro 
et al. 1999). Damping of these remnant short-trailing 
FMDWs inside the "ring" would induce a net accretion 
of gas and magnetic flux to fuel nucleus activities. For 
NGC 1097, in addition to the starburst "ring", circum- 
nuclear tight-winding trailing spiral arms of FMDWs are 
detected (Fig. 1) down to the immediate environs of the 
nucleus (r < 1.5"). We identify this spiral structure in- 
side the "ring" (Fig. 1) as remnant short-trailing FMDWs 
and suspect that the damping of these FMDWs plays a 
key role in fueling the central SMBH on sub-kpc scales. It 
is highly desirable that dedicated multiwavelength obser- 
vations of circumnuclear regions of spiral galaxies would 
provide more examples to enrich the scenario advanced 
here. While already indicative, it is anticipated that re- 
fined maps of polarized radio-continuum emissions would 
ultimately establish a more specific correspondence of spi- 
ral structures of gas density and magnetic field around 
and within the circumnuclear starburst "ring" , analogous 
to those manifest in galactic spiral structures on larger 
scales. We expect the total and polarized radio-continuum 
spiral arms to closely follow the dust lanes along the inner 
edges of circumnuclear optical spiral arms of NGC 1097. 



Fig. 1. — A wavelet transformed and reconstructed WFPC2/HST 
(F555W filter) image in false color of the central spiral arms around 
and within the circumnuclear starburst "ring" of the barred spiral 
galaxy NGC 1097 showing an angular scale lower than 2" . 

The large pitch angle of ~ 50° of optical arms (Rickard 
1975) and magnetic field (Beck et al. 1999) outside the 
"ring" is interpreted in terms of long-trailing FMDWs sus- 
tained by cf) E at the mILR. As a test, it may be possible to 
detect a weaker short-trailing FMDW outside the "ring" 
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